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Abstract. Since 2004, the MAGIC y-ray telescope has newly discovered 6 TeV blazars. The total set of 13 MAGIC-detected 
active galactic nuclei includes well-studied objects at other wavelengths like Markarian 501 and the giant radio galaxy M 87, 
but also the distant the flat-spectrum radio quasar 3C 279, and the newly discovered TeV y-ray emitter S5 07 16+7 1 . In addition, 
also long-term and multi-wavelength studies on well-known TeV blazars and systematic searches for new TeV blazars have 
been carried out. Here we report selected highlights from recent MAGIC observations of extragalactic TeV y-ray sources, 
emphasizing the new physics insights MAGIC was able to contribute. 
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INTRODUCTION 

The study of very high energy (VHE, E > 70 GeV) y-ray 
emission from active galactic nuclei (AGN) is one of the 
major goals of ground-based y-ray astronomy. The sen- 
sitivity of the current imaging air Cerenkov telescopes 
(IACT) enables studies of the blazar phenomenon, and 
in particular advances in understanding both the origin 
of the VHE y-rays as well as the relations between pho- 
tons of different energies (from radio to VHE). 

Except for the radio galaxy M87, all 23 currently 
known VHE y-ray emitting AGNs [1] are high-frequency 
peaked BL Lac objects [2], 1 a subclass of blazars char- 
acterized by a low high-energy luminosity when com- 
pared with quasars, and a synchrotron peak in the X- 
ray band. Their Spectral Energy Distribution (SED) is 
characterized by a second peak at very high y-ray ener- 
gies. In synchrotron-self-Compton (SSC) models it is as- 
sumed that the observed y-ray peak is due to the inverse- 
Compton (IC) emission from the accelerated electrons 
up-scattering previously produced synchrotron photons 
to high energies [3]. In hadronic models, instead, inter- 
actions of a highly relativistic jet outflow with ambient 
matter, proton-induced cascades, or synchrotron radia- 
tion off protons, are responsible for the high-energy pho- 
tons. Another defining property of blazars is the high 
variability of their emission ranging from radio to y-rays. 
For VHE y-ray blazars, correlations between X-ray and 
y-ray emission have been found on time scales ranging 
from w 10 minutes to days and months (see, e.g., Fos- 
sati et al. [4], Lenain [5]), although the correlations have 
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proven to be rather complicated [6] . 

Here we present selected results on multi-wavelength 
campaigns with MAGIC participation and for MAGIC 
blazar observations of Mkn 501, M87 (February 2008), 
1ES 1011+496, S5 0716+71, and 3C279. 



THE MAGIC TELESCOPE 

MAGIC [7], located on the Canary Island of La Palma 
(2200 m a.s.L), is currently the largest (17-m) single- 
dish IACT. Its energy range spans from 50-60 GeV (trig- 
ger threshold at small zenith angles) up to tens of TeV. 
MAGIC has a sensitivity of ~ 1 .6% of the Crab Nebula 
flux in 50 observing hours. Its energy resolution is about 
30% above 100 GeV and about 25% from 200 GeV on- 
wards. The MAGIC standard analysis chain is described, 
e.g., in Albert et al. [8]. Observations during moderate 
moonshine enable a substantially extended duty cycle 
[9], which is particularly important for blazar observa- 
tions. Parallel optical /?-band observations are performed 
by the Tuorla Blazar Monitoring Program and its KVA 
35-cm telescope. A second MAGIC telescope is being 
commissioned. By stereoscopic observations, the sensi- 
tivity of the two-telescope observatory will be further in- 
creased considerably. 

STRONG FLARING OF MESSIER 87 IN 
FEBRUARY 2008 

M87 is the only non-blazar radio galaxy known to 
emit VHE y-rays, and one of the best-studied extra- 
galactic black-hole systems. To enable long-term stud- 
ies and assess the variability timescales and the loca- 
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FIGURE 1. The night-by-night light curve for M 87 as measured from 2008 January 30 (MJD 54495) to 2008 February 1 1 (MJD 
54507). Left: The flux in the energy bin 150 — 350 GeV, being consistent with a constant emission. Right: Integral flux above 
350 GeV; flux variations are apparent on variability timescales down to 1 day. The inset shows the light curve above 350 GeV in a 
40 min time binning for the night with the highest flux (2008 February 1). The vertical arrows represent flux upper limits (95% c.l.) 
for the nights with negative excesses. 



tion of the VHE emission in M87, the H.E.S.S., VER- 
ITAS, and MAGIC collaborations established a regular, 
shared monitoring of M 87 and agreed on mutual alerts in 
case of a significant detection [10]. Results from the en- 
tire monitoring campaign will appear in a separate paper. 
During the MAGIC observations, a strong signal of 8(7 
significance was found on 2008 February 1, triggering 
the other IACTs as well as Swift observations. For the 
first time, we assessed the energy spectrum below 250 
GeV [11, 12], which can be described by a power law 
with a spectral index of F = 2.30 ± 0.1 l sta t ± 0.20 syst . 
This relatively hard VHE spectrum is unique among 
the extragalactic VHE y-ray sources, which show either 
curved or softer spectra. We did not measure a high- 
energy cut-off, but found a marginal spectral hardening, 
which may be interpreted as a similarity to other blazars 
detected at VHE, where such hardening has often been 
observed [1]. 

Our analysis revealed a variable (significance: 5.6 a) 
night-to-night y-ray flux above 350 GeV, while no vari- 
ability was found in the 150-350 GeV range (Fig. 1). We 
confirm the E > 730 GeV short-time variability of M 87 
reported by Aharonian et al. [13]. This fastest variabil- 
ity At observed so far in TeV y-rays in M 87 is on the 
order of or even below one day, restricting the emission 
region to a size of R < At c8 = 2.6 x 10 15 cm = 2.65 
Schwarzschild radii (Doppler factor 8), and suggests the 
core of M 87 rather than the brightest known knot in the 
M 87 jet, HST-1, as the origin of the TeV y-rays. 

There exists no lower limit on the size of HST-1, 
and thus the flux variability cannot fully dismiss HST- 
1 as possible origin of the TeV flux. During the MAGIC 
observations, however, HST-1 was at a historically low 
X-ray flux level, whereas at the same time the luminosity 
of the M 87 core was at a historical maximum (D. Harris, 
priv. comm.). This strongly supports the core of M 87 as 



the VHE y-ray emission region. 

THE JULY-2005 FLARES OF MKN 501 

Mkn 501 (z = 0.034) is known to be a strong and vari- 
able VHE y-ray emitter. MAGIC observed Mkn 501 for 
24 nights during six weeks in summer 2005 [14]. In 
two of these (one with moon present), the recorded flux 
exceeded four times the Crab-nebula flux, and revealed 
rapid flux changes with doubling times as short as 3 min- 
utes or less. For the first time, short (« 20 min) VHE 
y-ray flares with a resolved time structure could be used 
for detailed studies of particle acceleration and cooling 
timescales. 

Interestingly the flares in the two nights behaved dif- 
ferently: While the 2005 June 30 flare is only visible in 
250GeV-1.2TeV, the 2005 July 9 flare is apparent in all 
energy bands (120 GeV to beyond 1.2 TeV). Addition- 
ally, under some assumptions a time delay between the 
flare peak in different energy bins (E < 0.25 TeV and 
E > 1.2 TeV) in the latter flare was determined to be 
4 ± 1 minutes [14]. While this comparatively high value 
may be explained by the crude binned analysis, the indi- 
cated error is clearly underestimated. The reanalysis of 
these data in [15] resulted in a much-improved estimate 
for the time-energy relation. At a zero-delay probability 
of P = 0.026, a marginal time delay of T, = (0.030 ± 
0.012) sGeV~' towards higher energies was found using 
two independent analyses, both exploiting the full statis- 
tical power of the dataset (see [15, 16] for details). 

Several explanations for this delay have been consid- 
ered up to now: 

1. Particles inside the emission region moving with 
constant Doppler factor need some time to be ac- 
celerated to energies that enable them to produce y 



rays with specific energies [14]. 

2. The y-ray emission has been captured in the initial 
acceleration phase of the relativistic blob in the jet, 
which at any point in time radiates up to highest y- 
ray energies possible [17, 18]. 

3. An one-zone SSC model, which invokes a brief 
episode of increased particle injection at low ener- 
gies [19]. 

4. An energy-dependent speed of photons in vacuum 
[20], as predicted in some models of quantum grav- 
ity [21]. When assuming a simultaneous emission 
of the y-rays (of different energies) at the source, a 
lower limit of M QG > 0.21 x 10 18 GeV (95% c.l.) 
can be established [15], which increases further if 
any delay towards higher energies in the source is 
present. 

BLAZARS DETECTED DURING 
OPTICAL OUTBURSTS 

MAGIC has been performing target of opportunity ob- 
servations upon high optical states of known or poten- 
tial VHE blazars. Up to now, this strategy has proven 
very successful with the detection of Mkn 180, 1ES 
1011+496, and recently S5 0716+71 (see Lindfors & 
Mazin [26] and references therein). The 18.7-h observa- 
tion of 1ES 1011+496 was triggered by an optical out- 
burst in March 2007, resulting in a 6.2-a detection at 
^>200GeV = (1.58 ±0.32) x lO^cm-V 1 [27]. An in- 
dication for an optical-VHE correlation is given, in that 
in spring 2007 the VHE y-ray flux is >40% higher than in 
spring 2006, where MAGIC observed the blazar as part 
of a systematic search for VHE emission from a sample 
of X-ray bright {F\ k e v > 2 jJ, Jy) northern HBLs at mod- 
erate (z < 0.3) redshifts [28]. 

In April 2008, KVA observed a high optical state of 
the blazar S5 0716+71, triggering MAGIC observation, 
which resulted in a the detection of a strong 6.8-c signal, 
corresponding to a flux of F>4ooGeV ~ 10 _11 cm _2 s _1 . 
The MAGIC observation time was 2.6 h. The source was 
also in a high X-ray state [29]. 

The determination of the before-unknown redshifts of 
1ES 101 1+496 (z = 0.21, [27]) and S5 0716+71 (z = 
0.31, [30]) makes these objects the third-most and 
second-most distant TeV blazars after 3C279, respec- 
tively. 

MULTI- WAVELENGTH CAMPAIGNS 

For an advanced understanding of blazars, coordinated 
simultaneous multi-frequency observations are essential, 
yielding SEDs spanning over 15 orders of magnitude. 



MAGIC participated in a number of multiwavelength- 
campaigns carried out on five Northern-hemisphere 
blazars, which involved the X-ray instruments Suzaku 
and Swift, the y-ray telescope H.E.S.S. and MAGIC, and 
other optical and radio telescopes. 

• Mkn 421 was detected in a campaign April 2006 in 
all wavelengths (publication under preparation), 

• The observations of Mkn 501 in July 2006 revealed 
the lowest X-ray and VHE state ever observed. No 
variability in VHE y-rays was found, while an over- 
all increase of about 50% during one day was seen 
in X-rays. A one-zone SSC model describes this 
quiescent state of Mkn 501 well [22]. 

• Also campaigns on 1ES 1218+304 and 
1H 1426+428 have been carried out, during both 
of which significant X-ray variability has been 
observed. The VHE data are being analyzed. 

■ The VHE emission of PG 1553+113 showed no 
variability during the first multi-wavelength cam- 
paign on this blazar in July 2006 [23, 24]. 

■ 1ES 1959+650 showed VHE data among the low- 
est flux state observed from this object, while at 
the same time a relatively high optical and X- 
ray flux (both Swift/Suzaku) was found [25]. The 
SED could be modeled using similar parameters 
as needed for the SED measured in 2002, with a 
slightly more compact source and a slightly lower 
magnetic field. 

Further campaigns have been organized. 

DETECTION OF THE FLAT-SPECTRUM 
RADIO QUASAR 3C 279 

Observations of 3C 279, the brightest EGRET AGN 
[31], during the WEBT multi-wavelength campaign [32] 
revealed a 5. 77- a post-trial detection on 2006 Febru- 
ary 23 supported by a marginal signal on the preced- 
ing night [33, 34]. The overall probability for a zero-flux 
lightcurve can be rejected on the 5. 04- a level. Simulta- 
neous optical /?-band observations, by the Tuorla Obser- 
vatory Blazar Monitoring Program revealed that during 
the MAGIC observations, the y-ray source was in a gen- 
erally high optical state, a factor of 2 above the long-term 
baseline flux, but with no indication of short time-scale 
variability at visible wavelengths. The observed VHE 
spectrum can be described by a power law with a differ- 
ential photon spectral index of a = 4. 1 ± 0.7 stat ± 0.2 syst 
between 75 and 500 GeV (Fig. 2). The measured inte- 
grated flux above 100 GeV on 23 February is (5.15 ± 
0.82 stat ± 1.5 syst ) x 10 -10 photons cm" 2 s" 1 . 

Located at z = 0.536, VHE observations of such dis- 
tant sources were until recently impossible due to the ex- 
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FIGURE 2. Differential energy spectrum of 3C 279. The 
grey area includes the combined statistical (la) and systematic 
errors. The dotted line shows compatibility of the measured 
spectrum with a power law of photon index a =4.1. The tri- 
angles are measurements corrected on the basis of two models 
for the EBL density (see text). 



pected strong attenuation of y rays by the extragalactic 
background light (EBL; see [35] for a review), which in- 
fluences the observed spectrum and flux, resulting in an 
exponential decrease with energy and a cutoff in the 7- 
ray spectrum. The observed VHE spectrum is sensitive 
to the EBL between 0.2 — 2/im. The reconstructed intrin- 
sic spectrum is difficult to reconcile with models predict- 
ing high EBL densities (e.g., the fast-evolution model of 
Stecker et al. [36]), while low-level models, e.g. [37, 38], 
are still viable. Assuming a maximum intrinsic photon 
index of a* = 1.5, an upper EBL limit is inferred, leav- 
ing a small allowed region for the EBL. 

The results support, at higher redshift, the conclusion 
drawn from earlier measurements [38] that the observa- 
tions of the Hubble Space Telescope and Spitzer correctly 
estimate most of the light sources in the Universe. The 
derived limits are consistent with the EBL evolution cor- 
responding to a maximum star-formation rate at z > 1 , as 
suggested by [37] and similar models. 
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